Globally coupled bistable elements as a model of group decision making 



oo 



> : 
\o . 
o ■ 
o ■ 
in : 
o . 

On ■ 
On ; 

o 
in: 

& 
c3 



Damian H. Zanette 
Consejo Nacional de Investigaciones Cientificas y Tecnicas 
Centro Atomico Bariloche and Instituto Balseiro 
8400 Bariloche, Rio Negro, Argentina 

A simple mathematical model is proposed to study the effect of the average trend of a population 
on the opinion of each individual, when a group decision has to be made by voting. It is shown 
that if such effect is strong enough a transition to coherent behavior occurs, in which the whole 
population converges to a single opinion. This transition has the character of a first-order critical 
phenomenon. 



Group decision making is a complex social process in 
which the inherent factors that determine the position of 
each individual -such as previous experience, prospective 
benefits, current personal circumstances, and character- 
interact in a nontrivial way with the average trend, to 
which individuals are exposed through communication 
between them. Group decision results from this interac- 
tion as an emerging property of their collective behavior. 

Consider as a specific case an ensemble of individuals 
that have to choose, at a given time in the future and 
by individual voting, and option among a prescribed set 
of instances. After vote counting, the decision is sim- 
ply taken for the most voted option. This decision - 
which, once votes have been emitted, is straightforwardly 
defined- results however from the complex collective pro- 
cess that builds up the opinion of each individual [Q. 
During a certain period previous to the voting act, in 
fact, the opinion of a given individual evolves due to the 
modulation that the knowledge of another's position im- 
poses on the own tendency. In an efficiently communicat- 
ing ensemble, like in any modern population, individuals 
are continuously exposed to the average opinion of the 
ensemble -for instance, through poll results, published 
by mass communication media- and are expected to be 
more or less strongly influenced by this collective element 

This paper is aimed at exploring, in the frame of a sim- 
ple mathematical model, how the effects of the personal 
trend and of the average opinion in defining the individ- 
ual vote combine with each other to lead the group to its 
collective decision. For the sake of concreteness, suppose 
that the population has to choose by voting among two 
candidates, C + and C~ . In the model, the time evolu- 
tion of the opinion of the i-th individual is described by 
a variable Xj(i), with Xj G [ — 1,1] for all i and t. Large 
values of Xj, |xi| ~ 1, are to be associated with a firm 
decision to vote for one of the two candidates, (C + for 
x, > and C~ for x, < 0, say), whereas small values of 
x, correspond to a looser opinion. In any case, when the 
voting act takes place the individual opinion is quenched 
and the decision is made according to the sign of Xi at 
that time. In practice, it is supposed that the typical evo- 
lution times for the individual opinion are shorter than 
the time elapsed up to the voting act, so that one will 



focus the attention on the long-time asymptotics of the 
model. 

The average opinion, which is expected to play a rel- 
evant role in the definition of the individual decision, is 
here characterized by the arithmetic mean value 



(1) 



where N is the size of the population. This mean value 
is a measure of how much defined is the global trend to- 
wards on the two candidates. In fact, the sign of x at the 
time of the voting act determines the chosen candidate. 

To stress the effect of the average opinion on the in- 
dividual vote it is assumed that, in the absence of such 
effect, each individual would simply reinforce his or her 
initial personal opinion as time elapses. This means, in 
particular, that a given individual would not change his 
or her original preference for one of the candidates. This 
behavior is well represented by the following dynamical 
equation for Xj(t): 



dxi 
~dt 



(i = l,2,...,N). 
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In fact, the solution to this equation approaches the 
asymptotic value Xj(oo) = +1 or Xj(oo) = —1 depend- 
ing on the initial value Xi(0) being positive or negative, 
respectively. Moreover, Xj(i) does not change its sign 
during the whole evolution. If Xi(0) = 0, then xt(t) = 
for all t, but this stationary state is unstable. These facts 
can readily be verified from the explicit solution to Eq. 
(0), which reads, 



Xi(t) 



Xi(0) 



v / x l (0) 2 - \x0? ~ 1] exp(-2i) 



(3) 



From a dynamical viewpoint, Eq. (Q) implies that each 
individual behaves as a bistable element, its asymptotic 
state being fixed by the initial condition. In physics, this 
kind of model has been used to study spin systems 
the soft-spin approximation) |j| and neural networks 

The effect of the average opinion on the evolution of 
the individual trend is described by modifying Eq. (j~ 
the following way: 
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— i = a* - x f + ki(x - Xi) (i = l,2,...,N), (4) 



where x(t) has been defined in ([!]) and ki is a constant 
that, as discussed in the following, measures the influ- 
ence of the average opinion on the i-th individual. For 
hi > the new terms drive Xi(t) towards the average 
x(t). In fact, for large values of ki and slowly varying x, 
the individual variable Xi would exponentially approach 
the average. The new terms thus represent, for positive 
ki, a trend of the i-th individual to follow the average 
opinion, which can either reinforce or compete with his 
or her individual position. Negative values of ki would 
correspond to individuals who tend to take a position 
opposite to the average. 

In a physical context, the new terms represent an "in- 
teraction" between individuals. From a mathematical 
viepoint, in fact, those terms couple the set of equations 
through the average x, which depends on the whole 
set of Xi (i — 1, . . . , N). This coupling makes it impossi- 
ble to give the exact solution to the model equations (|4|), 
and the system has to be treated numerically. In partic- 
ular, note that it is not possible to derive an autonomous 
equation for the evolution of the average x(t). 

The case where the coupling constant ki is positive 
and the same for all individuals, fc; = k > for all i, 
is considered first. Obviously, for k — the uncoupled 
ensemble -whose behavior has been discussed above- is 
recovered. For k = 1, Eq. (||) reduces to 



Let Tjj = Xi — Xj be the difference between the states of 

... l~K 
any two individuals. It can be shown from Eq. (|3|) that, 

if — 1 < Xi, Xj < 1, Tij tends to zero as time elapses ||. 

In other words, for k = 1 the model predicts that all the 

individuals will have the same opinion at sufficiently long 

times. 

Figure 1 shows the evolution of Xi in a population of 
1000 individuals for k = and k = 1. For the sake of 
clarity, only 100 variables are displayed. Initially, the in- 
dividual states are randomly distributed in [—1,1]. As 
expected, for k — the states are soon divided into two 
clusters, according to their signs. For k = 1, instead, 
all states are attracted to a single cluster. Since for this 
value of the coupling constant the average opinion is al- 
ready dominant and all the individuals behave in a coher- 
ent way, it can be predicted that for k > 1 the dynamics 
of the ensemble is qualitatively the same. This is indeed 
verified from numerical results. On the other hand, for 
< k < 1 a transition is expected to occur between the 
two qualitative different behaviors observed at k = and 
k = 1. This transition is characterized in the following. 

According to numerical calculations, for sufficiently 
small values of the coupling constant the collective be- 
havior qualitatively reproduces the evolution of the un- 
coupled ensemble (k = 0). In fact, if the initial distribu- 
tion of Xi is uniform over [—1,1] the population becomes 
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FIG. 1. Evolution of Xi for a population of 1000 individ- 
uals with k — and k = 1. For clarity, only 100 states are 
displayed. The initial distribution is uniform over [—1, 1]. 

divided into two groups -as when, for k = 0, both signs 
are initially present. If, instead, one of the signs is much 
more abundant than the other, all the variables may ul- 
timately converge to one of the extreme values Xi = ±1 
-as when, for k = only one sign is initially present. 
The interacting ensemble is therefore "bistable," in the 
sense that two qualitatively different asymptotic states 
are observed depending on the initial condition: either 
all individuals behave coherently, or they become divided 
into two groups. On the other hand, as stated above, for 
larger values of k only coherent behavior is observed. Fig- 
ure 2 illustrates these behaviors for intermediate values 
of k. 

The transition between bistable and coherent behavior 
is due to a stability change in the possible asymptotic 
states of the coupled ensemble. Suppose that, as the 
system evolves, the N individuals become divided into 
two groups. One of them, with pN individuals (0 < 
p < 1) approaches the asymptotic state X\, whereas the 
other, with (1 — p)N individuals, approaches X^- It has 
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FIG. 2. Evolution of Xi for a population of 1000 individuals 
with k — 0.4 and k = 0.8. For clarity, only 100 states are 
displayed. The initial distribution is uniform over [—1,1]. 

to be stressed that the value of p depends in a nontrivial 
way on the initial condition, and cannot be analytically 
determined a priori. According to Eq. (|j) the following 
identities should hold for t — > oo and N — > oo: 

= (1 - k)Xi + k[ P X 1 + (1 - p)X 2 ] - XI 

(6) 

= (1 - k)X 2 + k[ P X 1 + (1 - p)X 2 ] - xl 

These equations include also the case of coherent behav- 
ior, if one puts X\ — X 2 with any value of p. Their 
solutions constitute the set of stationary states for the 
system, whose stability can be studied by means of stan- 
dard linearization around equilibria. 

Equations (g) can be reduced to a 9th-degree polyno- 
mial equation for either X\ or X 2l and have therefore 
nine solutions -which in general are complex numbers. 
The trivial solution X\ = X 2 — is unstable. The re- 
maining eight solutions can be grouped into symmetrical 
pairs, (Xi,X 2 ) and {—Xi, —X 2 ), both with the same sta- 
bility properties. It is therefore enough to analyze, for 
instance, the four solutions with X\ > 0. (i) The first 
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FIG. 3. The equilibrium states X\ and X2 as a function 
of k for solutions (ii), (iii), and (iv) (see text), at p = 0.52. 
For the solutions displayed, Xi is always positive and X2 is 
always negative. 

one, Xi — X 2 = 1, is stable for all k and corresponds to 
the asymptotic state of coherent evolution, (ii) The sec- 
ond solution is real and unstable for all k. It approaches 
the unstable solution (0,-1) for k — » and the triv- 
ial solution for k — > 1. (iii) Another unstable solution 
approaches the unstable state (0, 1) as k — » 0. (iv) Fi- 
nally, there is a stable solution that approaches (1,-1) 
as k — > 0. This solution corresponds to the case where 
the individuals have become divided into two groups. 

Figure 3 shows the numerical results for X\ and X 2 
as a function of k, for p = 0.52. Solid lines indicate 
stable solutions whereas dashed lines stand for unstable 
solutions. As the coupling constant grows, there is a 
critical value k c at which the two solutions (iii) and (iv) 
collide and become complex. At this critical value, the 
solution where the population is divided into two groups 
dissapears. The value of k c is related to p by the following 
polynomial equation: 

4fc c - 18k 2 c {p - p 2 ) + 21k\(p 2 - 2p 3 +p 4 ) = l. (7) 

Thus, for a given value of p -which is determined by the 
initial condition- and k < k c , two qualitatively different 
behaviors can occur. Either X\ = X 2 = ±1, and the sys- 
tem evolves coherently, or X 2 ^ X%, and the individuals 
are divided into two groups. For k > k c , instead, only 
coherent behavior is possible. Thus, for sufficiently large 
k, the opinion of the whole population approaches the 
same state. Figure 4 shows a phase diagram k versus p, 
where the boundary between the zones of bistability and 
coherence given by Eq. (]?]) is shown. 

The transition between bistability and coherence can 
be characterized by a single order parameter introducing, 
for instance, the mean difference 6 between the states of 
any pair of individuals, 

5= l -p{l-p)\X 1 -X 2 l (8) 
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FIG. 4. Phase diagram in the (p, fc) plane, displaying the 
zones of coherent behavior and bistability. 

which has been plotted in Fig. 5 as a function of k for 
p = 0.52. The dependence of 5 on the coupling constant 
suggests classifying the transition as a first-order critical 
phenomenon. 

Naturally, the assumption that all the individuals in 
the population are equally influenced by the average 
opinion -i.e. that all the individuals have the same cou- 
pling constant k— is not a very realistic one. Rather, it 
is to be expected that the coupling constants are dis- 
tributed within a certain interval, with some individuals 
being more influenced by another's opinion than other. 
It could moreover be supposed that some individuals are 
negatively affected by the mean trend, tending to make 
their own opinion diverging from the average. This case 
would correspond to fcj < 0. 

In this case of inhomogeneous behavior it can be easily 
shown from Eq. (^) that the asymptotic state of Xi(t) 
depends on the value of ki. This relation is implicitly 
given by the equation 

3 

k = (9) 

Numerical simulations show however that, as far as the 
distribution of coupling constants fc, is moderately nar- 
row, the qualitative collective behavior is the same as for 
uniform k. Mathematically, it is not an easy task to char- 
acterize the situation in which this behavior breaks down 
as the values of ki become more and more scattered. Is 
is nevertheless expected that coherent evolution can be 
destroyed if the coupling constants are sufficiently differ- 
ent from each other, including in particular some nega- 
tive values. A sufficient condition for coherence to fail is 
in fact that a single individual has a coupling constant 
ki < —2. In this situation, however, it is this only indi- 
vidual who fails to behave coherently, thus not affecting 
the result of the collective decision making. 

In summary, it has been here shown within a sim- 
ple mathematical model that, under a sufficiently strong 
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FIG. 5. The order parameter 5 as a function of k, for 
p — 0.52. At k c ~ 0.592, S vanishes suddenly. 

influence of the average trend of a population on the 
opinion of each individual, the population behaves co- 
herently and votes converge towards a single candidate. 
For weaker coupling between individuals, instead, votes 
are more evenly divided between the two candidates. 
The transition between both behaviors is abrupt and, 
in fact, has the character of a critical phenomenon. This 
is qualitatively similar to the ferromagnetic phase transi- 
tion observed in spin systems -though in this physical 
phenomenon the phase transition is of the second or- 
der ||. Indeed, beyond the critical point the state of 
all the elements in the system coincide, even in spite of 
the initial condition corresponding to a uniform distribu- 
tion of states. The coupling mechanism is in fact able 
to break the initial macroscopic homogeneity, enhancing 
microscopic fluctuations. It can be interesting to further 
analyze this model, including for instance local communi- 
cation ways between individuals as well as noise, that can 
perturb in a nontrivial way the properties of the quoted 
transition 

The present results could encourage unfair, unscrupu- 
lous candidates to manipulate poll results published in 
mass media -if they have the power to do so- in their 
own benefit. 
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